home *** CD-ROM | disk | FTP | other *** search
/ Enter 2004 January / enter-2004-01.iso / files / maxima-5.9.0.exe / {app} / share / maxima / 5.9.0 / src / numerical / slatec / zbunk.lisp < prev    next >
Encoding:
Text File  |  2003-02-09  |  1.5 KB  |  44 lines
  1. ;;; Compiled by f2cl version 2.0 beta 2002-05-06
  2. ;;; 
  3. ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
  4. ;;;           (:coerce-assigns :as-needed) (:array-type ':simple-array)
  5. ;;;           (:array-slicing nil) (:declare-common nil)
  6. ;;;           (:float-format double-float))
  7.  
  8. (in-package "SLATEC")
  9.  
  10.  
  11. (defun zbunk (zr zi fnu kode mr n yr yi nz tol elim alim)
  12.   (declare (type (simple-array double-float (*)) yi yr)
  13.            (type f2cl-lib:integer4 nz n mr kode)
  14.            (type double-float alim elim tol fnu zi zr))
  15.   (prog ((ax 0.0) (ay 0.0))
  16.     (declare (type double-float ay ax))
  17.     (setf nz 0)
  18.     (setf ax (* (abs zr) 1.7321))
  19.     (setf ay (coerce (abs zi) 'double-float))
  20.     (if (> ay ax) (go label10))
  21.     (multiple-value-bind
  22.         (var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-7 var-8 var-9 var-10
  23.          var-11)
  24.         (zunk1 zr zi fnu kode mr n yr yi nz tol elim alim)
  25.       (declare
  26.        (ignore var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-7 var-9 var-10
  27.         var-11))
  28.       (setf nz var-8))
  29.     (go label20)
  30.    label10
  31.     (multiple-value-bind
  32.         (var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-7 var-8 var-9 var-10
  33.          var-11)
  34.         (zunk2 zr zi fnu kode mr n yr yi nz tol elim alim)
  35.       (declare
  36.        (ignore var-0 var-1 var-2 var-3 var-4 var-5 var-6 var-7 var-9 var-10
  37.         var-11))
  38.       (setf nz var-8))
  39.    label20
  40.     (go end_label)
  41.    end_label
  42.     (return (values nil nil nil nil nil nil nil nil nz nil nil nil))))
  43.  
  44.